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1.0  SUMMARY 

Van  Driest's  theory,  which  provides  a  continuous  velocity  and  shear 
distribution  for  turbulent  flow  near  a  nonporous  wall,  is  extended  to 
turbulent  flow  near  a  porous  wall.  Bie  new,  modified  theory  enables  the 
theoretical  calculation  of  velocity  profiles  to  be  performed  for  a  wider 
range  of  mass-transfer  rates,  and  it  gives  good  agreement  with  experimental 
data. 
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4.0 


PRINCIPAL  NOTATION 
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5.0  INTRODUCTION 


In  reference  1,  Van  Driest  introduced  a  very  useful  modification  of 
Prandtl's  mixing  length  theory,  which  provided  a  continuous  velocity  and 
shear  distribution  for  flat-plate  turbulent  flow  near  a  nonporous  wall. 

This  modification  also  formed  the  basis  for  the  theoretical  calculation  of 
the  velocity  profiles  and  has  been  used  quite  successfully  by  several 
investigators,  for  example  by  Cebeci  and  Smith  [2],  and  Patankar  and 
Spalding  [3], 

Hie  purpose  of  this  report  is  to  show  a  possible  method  of  extending 
Van  Driest's  modification  to  turbulent  flow  near  a  porous  wall.  The  new 
modification  enables  the  theoretical  calculation  of  velocity  profiles  to 
be  performed  for  a  wide  range  of  mass-transfer  rates  and  gives  good  agree¬ 
ment  with  experimental  data. 

It  is  important  to  note  that  the  present  modification  increases  the 
capability  of  the  Douglas  boundary-layer  method  and  allows  the  calculation 
of  turbulent  flows  with  mass  transfer,  which  becomes  quite  important  in 
boundary-layer  control  and  drag-reduction  problems,  as  well  as  in  problems 
where  it  is  necessary  to  protect  surfaces  exposed  to  high-temperature  gases. 


6.0  ANALYSIS 


6.1  Van  Driest 'a  Analysis 


Consider  an  infinite  flat  plate  undergoing  simple  harmonic  oscillation 
parallel  to  the  plate  in  an  infinite  fluid.  As  was  shown  by  Stokes,  the 
amplitude  of  the  motion  diminishes  from  the  wall,  as  a  consequence  of  the 
factor  exp(-y/A),  where  A  is  a  constant  depending  upon  the  frequency  of 
oscillation  of  the  plate  and  the  kinematic  viscosity  of  the  fluid.  Hence, 
when  the  plate  is  fixed  and  the  fluid  oscillates  relative  to  the  plate,  the 
factor  [1  —  exp(-y/A)]  must  be  applied  to  the  fluid  oscillation  to  obtain 
the  damping  effect  of  the  wall.  By  using  this  ariiunent,  Van  Driest  proposed 
that  the  mixing  length  l  in  Prandtl's  mixing-length  theory, 


■x 


t 


should  be  changed  to 


l  =  *y[l  -  exp(-y/A)] 


(6.1) 


(6.2) 


because  of  the  presence  of  the  wall,  since  l  -  <y  does  not  hold  in  the 
viscous  sublayer.  In  dimensionless  quantities  Eq.  (6.2)  can  also  be  written 


l 


*y[l  “  exp(-y+ /A+)] 


(6.3) 


Taking  k  =  0.4  and  A+  =  26,  Van  Driest  observed  that  calculated  velocity 
profiles  near  the  wall  correlated  very  well  with  the  experimental  data.  As 
a  result,  he  proposed  that  the  constant  A  should  be 


A 


(6.4) 
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6.2  lYeocnt  Aimlyttia 


lhc  expression  given  by  E<j.  (C. ^ )  w no  obtained  l>r  a  Mat-plate  M  v  with 
no  mass  trarsfor.  As  it  stands,  It  cannot  b©  used  !<»r  Mows  with  pressure 
gradienta.  Thia  ia  quite  obvious,  since  for  a  flow  with  an  adverse  pressure 
gradient,  i  -*  0,  which  Introducen  a  discontinuity  it>  mixing- length  cxpres- 
aion  and  conaequently  in  velocity  profiles.  It  eon  be  extended  to  M-jws  with 
preaaure  gradienta  by  the  following  reasoning. 

Conalder  the  momentum  equation  for  tve-dlmensional  incompressible 
turbulent  flowa: 


In  the  laminar  aublayer  region  thia  equation  can  be  written  aa 


u 


-  i  <l£  ♦ 

p  dx 


1  M 
p  $y 


by  neglecting  the  Reynolds  ahear-atreaa  term  -pu'v',  and  by  putting 
i  =  p  du/dy.  If  the  expreaaion  (6.  **)  ia  written  aa 


A 


('•  7) 


In  other  words,  if  t  is  treated  as  a  variable  and  not  aa  a  constant  in  the 
y-direction*,  then  for  the  region  close  to  the  wall,  wo  can  write  Eq.  (6.6)  as 


T 


(o.  (;) 


providing  that  the  wall  is  nonporous  (v  *  0).  In  that  ca3e,  tnc  expression 
for  A  should  be 


A  =  ?.6\ 


(••  0 


♦In  reference  3,  Patankar  and  Spalding,  arguing  that  the  local  damping  is  more 
likely  to  be  affected  by  the  local  shear  stress  than  the  damping,  at  the  extreme 
edge  of  the  boundary  layer,  also  used  the  local  value  op  i  in  Eq.  (o.t). 
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which  means  that  the  pressure-gradient  effect  changes  the  damping  constant  A 
to 


A*  =  26 


1  +  3E-1 


dx  p(u*)2 


-1/2 


Equation  (6.10a)  can  also  be  written  as 

A+  =  26  [  1  -  pVl V2 


(6. 10a) 


(6.10b) 


+  + 
where  p  is  an  acceleration  parameter  defined  as  p  =  —  dp/dx  v/p(u*)' 

+  +  Jt  • 

and  y  is  a  Reynolds  number  defined  as  y  =  yu  /v. 


In  the  case  of  a  porous  wal]*  Eq.  (6.6)  can  be  approximated  as 


dr _ w  x  dp 

dy  v  dx 


(6.11) 


The  solution  of  Eq.  (6.11)  is 


t  = 


i  t;  +  v  exp ( ir y) 


(6.12) 


which  leads  to 


A=26v  f 


%  [exp  (^y)  -1] 


1  I  +  (u*)2  exp 


(?>)) 


-1/2 


(6.13) 


As  a  result, the  damping  constant  A  for  a  porous  wall  becomes 


( 


A+  =  26/<£_^ 


■ dx  ,  „v2 
(  pvw(u*) 


j^exp(^  y)  -  lj  +  exp  y)| 


-1/2 


(6.14a) 
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Equation  (t>.  l*»a)  can  also  be  written  as 


A 


♦ 


(«..!«*) 


It  should  be  noted  that  Eq.  (6.  IU)  reduces  to  the  fom  given  by  Eq.  ( 6.  1C )  Tor 

pressure  gradient  flows  with  no  muss  transfer  and  to  the  numerical  value 

for  flat-plate  flows  with  no  mass  transfer,  as  originally  proposed  by  Van  Driest. 

It  is  easy  to  see  from  Eq.  (6.10)  that  the  tern  in  d rackets  may  be  nega¬ 
tive  for  highly  accelerating  flows.  In  such  cases,  since  square  root  is 
involved,  this  will  lead  to  a  numerical  difficulty  that  can  easily  be  elimi¬ 
nated  by  taking  the  absolute  value  of  the  tern  in  brackets.  This  was  the 
procedure  foliated,  for  example,  in  reference  2.  On  the  other  hand,  if  we 
were  to  identify  A*  as  a  damping  constant  for  a  given  atrcamwlae  location, 
as  Van  Driest  did,  then  the  above  procedure  is  not  correct.  That  is,  A*  in 
either  Eq.  (6.10)  or  (6.11*)  should  not  be  a  function  of  y+.  Consequently, 
it  is  necessary  to  take  y*  as  a  constant.  Since  Eq.  (6.6)  is  valid  only  in 
the  sublayer  region.  It  will  be  assumed  that  the  sublayer  Reynolds  number, 
y+,  does  not  deviate  appreciably  from  the  unblcwn  flat-plate  sublayer  value, 
and  it  will  be  taken  as  11.8.  Equation  (6.1^b)  then  bccnr.cs 


(  ♦  r  + 

1  ♦  )  -1/2 

]-  ^  [exp  (11.6  v  )  -  1 

♦  exp  (11.8  v  ) > 

(6. l^c) 

1  Vw  1- 

J  ) 

Figure  1  shews  the  deviation  of  the  damping  constant  from  that  of  a  flat- 
plate  flow  with  no  mass  transfer,  that  is,  A+/26.  Figure  la  shows  a  plot  of 
Eq.  (6.10).  According  to  this  plot,  the  damping  constant  of  a  pressure-gradient 
flow  deviates  considerably  from  that  of  a  flat-plate  flew  and  becomes  quite 
large  as  p+  approaches  the  value  0.08.  This  means  that  at  large  values  of 
p+,  a  condition  corresponding  to  highly  accelerating  flew,  the  sublayer 
becomes  larger.  In  this  cas;,  the  skin-friction  values  will  be  lower  than 
those  obtained  by  A+  =  26,  corresponding  to  flat-plate  flew.  Figure  lb 
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shows  a  plot  of  &j,  (6.  ll*c)  for  /.oro-pressurc-gradlcnt  flow.  It  can  be  seen 
that  as  v*  increases  the  damping  constant  decreases.  In  this  case  the 
region  where  Prandtl's  unmodified  mixing  length  applies  moves  closer  to  the 
wall  and  results  in  smaller  sublayer  thicknesses  than  those  of  flat-plate 
flows  without  mass  transfer.  Figure  lc  Is  a  plot  of  Eq.  (6.11*)  that  shows 
the  effect  of  v*  on  the  damping  constant  for  several  p+  values.  The 
results  Indicate  that  the  damping  constant  Increases  appreciably  In  acceler¬ 
ating  flows  (positive  p*)  and  with  suction  increases  even  mere,  as  Is  to  be 
expected.  Figure  1J  shows  the  effect  of  p+  on  the  damping  constant  for 
several  v+-vnlucs.  As  Is  to  be  expected,  the  damping  constant  increases 
with  suction  and  positive  p  and  leads  to  larger  sublayer  thicknesses. 
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7.0  COMPARISON  WITH  EXPERIMENT 


Needless  to  say,  the  proposed  modification  of  Van  Driest's  modified 
mixing-length  expression  Is  empirical,  and,  like  most  expressions  used  in 
calculations  of  turbulent  flows,  It  must  be  checked  with  experiment.  This 
was  dono  by  comparing  the  calculated  A+-valuos  from  Eq.  (6. 14)  with  those 
obtained  experimentally  and  by  using  the  eddy-visoosity  method  described  In 
reference  2.  In  the  latter  case  the  boundary  layer  is  regarded  as  a  composite 
layer  consisting  of  Inner  and  outer  regions  for  which  separate  expressions  for 
eddy  viscosity  arc  used.  In  other  words,  the  solution  of  Eq.  (6.5)  is  obtained 
by  eliminating  the  Reynolds  shear-stress  term  by  the  following  eddy-viscosity 
formulation: 


€ 


(7.1) 


where  e  is  given  oy 


€  = 


:i  *  <o-l‘y>2[i-«p  (-Jr)]  || 

o-°-0l6&'eKl[1  +  5-5(^)6]  1 


0  S  y  5  yc  (7.2a) 


yc  a  y  i  5  (7.2b) 


and  the  matching  point  y  between  two  regions  is  obtained  from  the  continuity 

4. 

of  eddy  viscosity,  namely,  -  c  .  In  Eq.  (7.2a)  the  damping  constant  A 
is  given  by  Eq.  (6. lU). 

Figure  2  shows  the  variation  of  damping  constant  A+  with  blowing 
parameter  V  .  The  experimental  values  of  A+  were  obtained  from 

the  data  of  Simpson  ot  al  [4]  and  Kendall  [5],  and  were  reported  in  reference  6 
by  Bushnell  and  Beckwith.  The  figure  also  shows  the  curve  faired  to  the 
experimental  data  used  by  Bushnell  and  Beckwith,  together  with  the  results 
obtained  from  Eq.  (6.1.4)  for  zero  pressure  gradient.  The  skin-friction  values 
for  Eq.  (6. l4)  were  obtained  from  Simpson's  data  [4]  for  blowing  and  from 
Tennekes'  data  [7]  for  suction.  The  agreement  between  Eq.  (6.l4),  the  experi¬ 
mental  data,and  the  faired  curve  is  very  good  for  blowing  parameters  up  to  l4. 

For  larger  blowing  parameters,  the  calculated  A+  values  deviate  from  the 
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faired  curve  uoed  by  Bushnell  and  Beckwith,  but  seem  to  agree  reasonably  wt  LI 
with  experimental  data  except  for  one  value. 


Figures  3a,  3b,  and  3c  show  the  calculated  velocity  profiles  obtained 

by  the  eddy -viscosity  formulation  given  in  Eq.  (7.2)  for  blowing  rates  of 

v  /u  =  0.00386,  0.00784,  and  0.00950,  respectively,  for  the  experimental  data 
we 

of  Simpson  et  al  [4],  In  these  calculations  the  damping  constant  A  as  given 
by  Eq.  (6. 14)  was  used.  It  is  important  to  note  that  calculations  could  not 
be  mpde  with  the  unmodified  damping  constant  given  by  Eq.  (6.4).  Figures  4a 
and  4b  show  a  comparison  of  calculated  and  experimental  skin-friction  values 
obtained  for  two  blowing  rates,  \/ne  =  0.00386  and  0.00784,  respectively. 

In  general,  the  agreement  seems  satisfactory,  and  a  slight  dicrepancy  in  skin- 
friction  values  could  be  due  to  the  procedure  that  was  followed  in  making  com¬ 
parisons.  Ideally,  one  should  start  the  theoretical  calculations  by  inputting 
the  initial  velocity  profile  obtained  from  the  experimental  data.  However,  in 
the  present  comparisons,  an  effective  length  that  matched  the  momentum  thickness 
at  the  station  where  blowing  began  was  determined.  This  procedure  does  not 
necessarily  match  the  initial  skin-friction  coefficient  at  that  station.  How¬ 
ever,  as  the  calculations  continue  in  the  streamwise  direction,  the  initial 
discrepancy  in  the  local  skin  friction  decreases.  For  example,  when  the  com¬ 
parisons  were  made  for  the  blowing  ratio  of  v  /u  =  0.00950,  the  calculated 

^  ^  —4  -4 

and  experimental  skin-friction  values  were  5.5  x  10  and  3*4  x  10  , 

respectively,  at  the  station  where  the  matching  was  made  (R^  =  4.3  x  10^). 

Downstream,  at  sufficiently  higher  Reynolds  number,  the  discrepancy  became  quite 

small.  For  example,  at  R-  =  15  x  10^,  the  experimental  and  calculated  local 

**  -4  —4 

skin  friction  values  were  1  x  10  and  1.1  x  10  ,  respectively,  and  at 

3  -4  -4 

Rg  =  17  x  10  ,  they  were  1  x  10  and  O.98  x  10  ,  respectively. 

Comparisons  were  also  made  for  the  sucked  boundary  layers.  The  experi¬ 
mental  data  of  Tennekes  were  used  for  this  purpose.  Figure  5  shows  a  compari¬ 
son  of  calculated  and  experimental  velocity  profiles  for  a  suction  rate  of 
vw/ue  =  -0.00312.  Figures  6a  and  6b  show  comparisons  of  calculated  and 
experimental  values  of  momentum  thickness  and  local  skin  friction,  respectively. 
In  general,  the  agreement  is  quite  satisfactory. 


12 


For  the  test  cases  considered  in  this  report,  the  calculated  results 
indicate  that  the  proposed  modification  of  Van  Driest 's  expression  is  quite 
useful  and  gives  good  agreement  with  experiment.  Further  comparisons  on  the 
same  subject,  which  will  be  reported  in  reference  8,  also  show  the  same  good 
agreement  and  give  further  support  to  the  proposed  modification. 
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8.0  CONCLUDING  REMARKS 


The  results  show  that  the  proposed  modification  of  Van  Driest' s  theory 
seems  to  be  suitable  for  carrying  the  calculations  right  down  to  the  wall  to 
compute  flows  with  pressure  gradients  and  mass  transfer.  In  the  absence  of  a 
t  tter  empirical  Tit,  the  new  formulation  appears  to  be  quite  useful. 
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Figure  1b.  Deviation  of  damping  constant  of  a  porous  flat-plate  flow  from  that  of  a 
nonporous  flat-plate  flow. 
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Figure  2.  Comparison  of  calculated  and  experimental  damping  constants  for  a  flat-plate 
flow  with  mass  transfer. 
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Figure  3.  Comparison  of  calculated  and  experimental  velocity  profiles  for  the  blown 
boundary  layer  measured  by  Simpson  et  al  (4] ;  (a)  v^/uj  ■  0.00386, 

(b)  -  0.00784,  (c)  vw/ue  ■  0.00950. 
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Figure  4.  Comparison  of  calculated  and  experimental  akin  friction  values  for  the  blown 
boundary  layer  measured  by  Simpson  et  al  [4] ;  (a)  v^u.  -  0.00386,  and 
(b)  Vyy/Ug  -  0.00784. 
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Figure  5. 
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Comparison  of  calculated  and  experimental  velocity  profiles  for  the  sucked 
boundary  layer  measured  by  Tennekes  [7]  for  v^Ug  *  —0.00312. 
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Figure  6.  Comparison  of  calculated  and  experimental  (a)  momentum  thickness  values, 
and  (b)  skin-friction  values  for  the  sucked  boundary  layer  measured  by 
Tennekes  [7]. 
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